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Abstract 

Flavor changing interactions of the gluino allow the b quark to decay into the 

strange quark plus a gluino pair if the gluino is in the ultra low mass window 

below 1 GeV. In this case the enhancement of the nonleptonic b decay could 

explain the anomalous semileptonic branching ratio. 
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In the last few years it has been noted by many authors that a light gluino would help 
to explain several anomalies at the Z scale and other discrepancies between experiment and 
theory 0-^]. Also, the impact of the light gluino on the branching ratio b — > sj has been 
recently investigated in Surprisingly enough, the region of gluino mass below 0.7 GeV 
is poorly constrained by experiment |J. Here, assuming that the gluino is in this low mass 
window we propose that the decay b — > sgg might contribute considerably to the b total 
width thus reducing the theoretical prediction of the B semileptonic branching ratio. 

Currently, the experimental value for the semileptonic branching ratio is BRsl(B) = 

exp 

(10.43 ± 0.24)% ||, while the theoretical prediction gives a lower bound of 12.5% 0. The 
last number includes also perturbative QCD corrections. Non-perturbative QCD corrections 
are not expected to increase the inclusive non-leptonic widths of B mesons significantly. 
Therefore, the discrepancy between theory and experiment, which is at least 14%, still has 
to be explained. 

One solution proposed in |7|,[8|] could be a T(b — > ccs) enhancement due to larger than 
expected non-perturbative corrections in the b — > ccs channel. This would, however, enhance 
also the charm multiplicity to about 1.3 which would be more than 15% higher than the value 
from current experimental data. Unless the future measurements of the charm multiplicity 
lead to its expected value, it could be plausible that there are new contributions to charmless 
b decays, that were unaccounted for in the theoretical prediction of BRsl{B). Some authors 
H prefer rather conservative explanations of the B semileptonic branching ratio puzzle and 
suggest that the most natural solution lies within QCD itself. Nevertheless, they do not 
exclude the possibility of scenarios from beyond the standard model contributing to the 
solution of the problem. 

In the following we show that the b decay to the s quark and a light gluino pair can 
easily increase the nonleptonic branching ratio by 20% in certain regions of its parameter 
space. We also show that this process contributes to the total width of b by a considerably 
larger amount than the process b — > sg if the gluino is below 1 GeV. 

The process we are dealing with is a tree-level flavor-changing-neutral-current process 



with a down type squark in the intermediate state. The decay rate is calculated using the 
quark-squark-gluino Lagrangian M, given by 
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where p stands for the quark generation (in our case p = b or s) and % labels the squark states 
[i = bL, 6r, sl, sr, cIl, d,R). The X a are the eight generators of color SU(3). The matrices Tl 
and T# are (6 x 3) matrices given by 
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where U is the matrix that diagonalizes the down- type squark mass matrix squared, Mj. 



M I is written as 
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Adopting the notation of | 

m\ L l + MJ + cK^MlK Am M d 
Am M d m 2 m I + Mj 

in a basis where the 3x3 down-type quark mass matrix is diagonal. The matrices M u 
and Md are diagonal up- and down-type quark mass matrices respectively and K is the 
Cabibbo-Kobayashi-Maskawa mixing matrix. For simplicity we take tuqi = mofl and equal 
to the universal scalar mass mo- In order to simplify the process of analytic diagonalization 
of Mj we take the trilinear scalar coupling A equal to zero, which does not affect the result 
significantly. The c-parameter, which is responsible for flavor-violating interactions, plays 
an important role in our numerical estimates of the b —>■ sgg branching ratio. Some authors 
T0|] take c of order 0.01 or even lower, while others || suggest that c can be somewhat larger 



in magnitude. As regards the sign of the c-parameter, c < is preferred in the MSSM. In 
this paper, we treat c as a phenomenological parameter to be experimentally constrained. 

In the case A = 0, Mj is block- diagonal and only the upper-left block needs to be 
diagonalized. The upper-left block can be written in the form 



(4) 







( 





o\ 




(° 





0\ 




1 + 











+ c> 








e 






u 





bj 




u 


e 


1/ 



3 



where we have neglected the masses of the d and s quarks with respect to the mass of the 
b quark, and similarly m u and m c with respect to m t . The modified parameter d is equal 
to cm 2 /ml and b is m\jm\. Only the two leading terms were kept in the product K^M^K, 
namely those proportional to \K t b\ 2 and K^ s K t b, the first being taken equal to unity and the 
latter being denoted by e. The matrix that diagonalizes -^/J/ 3x3 ) is found to be 
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where / is a function of the variables b, c' and e defined by 
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f(b,c',e) = V(fc + C) +4eV 2 . 
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The complete diagonalizing matrix U is given by 

U(3X3) 0' 

/ 

where / is the (3 x 3) identity matrix. The evaluation of [/'Mj[/ gives a diagonal matrix 
with the squark masses squared on the diagonal. For the left handed mass-eigenstates one 
gets 
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m l L = m l [ 1 + 2 + 2 + 2 nb,C ' ,e) ) , 
while the right handed ones get masses 
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The matrix T L needed for calculation of the decay rate T(b 
following way 
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Note that the matrix UJ 3x3 -j reduces to the identity matrix in the limit c — > 0, as it should. 
The matrix Tr is trivially found to be 
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Having found the exact form of the matrices Tl and Tr, the calculation of T(b — > sgg) can 
be completed analytically. The invariant matrix element Ai consists of terms corresponding 
to the exchange of &£,, bn, sl and sr. We have neglected much smaller terms with di or cIr 
exchange. After performing the 3-body Lorentz invariant phase space integration, the decay 
rate becomes 

r(6 - sgg) = 2^1 (^) £ -1_ (t^ + r*W) ^ + iW) , (14) 

for i,j = bL,bn, sl, sr. The function I{x) is given by 

I(x) = l-8x 2 + 24x 4 lnx + 8x 6 -x 8 . (15) 

The overall multiplicative factor of 2 in Eq. ([TJ]) is due to the Majorana nature of the 
external gluinos. Using the "diagonal" character of Tr this can be futher reduced to 

^-^) = ^©5^r- W ri«. (16) 

The sum can be written in terms of the squark masses 

i /r 2 r /2 ( rn~ — m, 2 - 1 

\ 3 mwr L L L ~ p « ■ { } 

Using the expressions for the squark masses obtained as the eigenvalues of Mj we can write 
the result for T(b —>■ sgg) in the form 

27n \m b J 



x 



/ 

c 



l + 4 +c 4_ c 2|X*|24 



(18) 



4 

Both the terms m 2 /mQ and c 2 1 i<Q* | 2 ^| in the denominator can be neglected with respect to 

2 

c^|, if m is larger than 80 GeV. Note also that T(b — > sgg) cannot develop a pole because 



of the experimental lower limit on the masses of squarks in the intermediate state. For 
example, according to |2j we can require in the light gluino case that my > 60GeV. As will 
be discussed below, this imposes an additional constraint on the c parameter as a function 

of TUq. 

It is convenient to define the ratio 

T(b — > sgg) 

s ' 9 ' 9 = T{b^cud) + T{b^cusy (19) 
where the denominator is given by 

Y(b — > cud) + Y(b — > cus) 

2,G 2 F ml\K cb \ 2 (ml \ 

The expression for the phase space factor Jo and the values of leading-log anomalous di- 
mension enhancement r\ and next-to-leading corrections enhancement J can be found in the 
literature (e.g. an d their product is of order 0(1). Combining equations (|T8|), ([19]) and 
(pO|) we get 

128 fa s \ 2 sm 4 (6 w )c 2 fm w m t \ 4 K* s 2 
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If the ratio R S gg gets as high as 20%, then the branching ratio of the non-standard 
model decay b — > sgg is more than 14%. This could completely account for the discrepancy 



exp 



between BRsl{B) and BRsl{B) . The current experimental data can be fit if R 



QCD 



sgg 



(22.68 ±0.52)%. 

In general, by requiring R sgg to have a certain value, one gets c as a function of the 
universal scalar mass m$. This function is plotted in Fig. 1 for several different values of 
R sgg - We have used as(m b ) = 0.18, a(m b ) = 1/133, sin 2 6^ = 0.232 and m t = 170GeV 
in all numerical calculations. Also, we have taken advantage of the equality \K ts \ — \K cb \. 
From the first figure it can be seen that the needed nonleptonic enhancement in b decays 
can be obtained using the contribution from the process b —> sgg for reasonable values of 
mo and c. The necessary values of c as a function of mo are intermediate between those 
considered by refs 0] and [pl| . 



6 



The lower bound on the b squark mass of 60 GeV mentioned above does not interfere 
with any of the curves plotted in Fig. 1. In fact, the mass of the b squark is certainly above 
75 GeV for m > 80GeV. The function of m Q obtained using equations flID"D and 

(pl|) is plotted in Fig. 2. This has an interesting implication for the problem of the b excess 
in Z decays. A possible explanation of the b anomaly could have been the Z decay into the 
b squark and the b anti-squark. For that to be possible, the b squarks would have to have 
a mass less than Mz/2. In order to explain the 6-anomaly using this process, one would 
need the ratio T(Z — > bb)/T{Z — > bb) to be about 2%. This is possible, however, only if 
m~ bL ~ 0A7Mz- This imposes the following constraint on the c-parameter through Eq. (|10"D 
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Unfortunatelly, this constraint is incompatible with the c dependence on itlq that we got 
from the analysis of the B semileptonic branching ratio. In addition such a light b squark 
would lead to unacceptably large contributions from Z — > bbg. Therefore, in the light gluino 
scenario one has to rule out the possibility of b squarks being lighter than Mz/2. Never- 
theless, there are other mechanisms, that can explain the b excess |J without contradiction 
with our current analysis. The current calculation is not sensitive to the gluino mass varying 
in the range of the low mass window (0 GeV to 0.7 GeV). 

It is interesting to compare the decay rate (|1^) to the decay rate of b — » sg, because 
processes like this could also account for the missing 14-20% in the hadronic branching ratio 
of the B H . We use the formula for Ysusyib — > sg) given in [JEJ, that corresponds to the 



processes with a squark and a gluino exchange within a loop and an external gluon attached 
either to the gluino line or the squark line. The decay rate is given by 
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with x = m^/m 2 and rh 2 = ~ fm| + mfj. The functions (7 and / are given in [12| but for 
the purpose of our comparison we need only their limits as x — > 0. These are 

1 1 

lim/ a (z) = -; ]imf b {x) = -. (24) 

The decay rate in the case of the gluino with a negligible mass is then equal to 

20736 m 8 

where the terms proportional to m 2 /ml and its higher powers were neglected. Dividing Eq. 
(HI) by Eq. ( |18D one gets 

Fsusvib -» gg) _ 289^5 
r(6 -> s^) 768 vr ' 1 j 

indicating that the contribution to the total b width from b — > s^^ decay is dominant over 
the one from 6 — > in the light gluino scenario. The mechanism of Jl2| for the non-leptonic 



enhancement is only consistent with a b squark above M^/2 if the gluino is heavier than 2 
GeV and m is less than 150 GeV. 

In conclusion we can say that, assuming the gluino is light, the decay b — ► sgg provides 
a plausible explanation of the gap between BR SL (B) and BR SL (B) 



cxp 



. The final 

QCD 

state gluinos in b decay could hadronize into the gluino-gluon or gluino-gluino bound states 
discussed in [f| or merely into intrinsic gluino components of normal hadrons. Since the 
values of the c-parameter and the universal scalar mass are not yet well determined, the 
branching ratio of b — > sgg may provide a useful constraint as experiments improve. 

We acknowledge useful and stimulating discussions with P. Cox. This work was supported 
in part by the U.S. Department of Energy under Grant No. DE-FG05-84ER40141. 
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FIGURES 

FIG. 1. The absolute value of c (c is assumed to be negative) is plotted as a function of the 
universal scalar mass mo for R S gg = 0.3 (dot-dashed line), 0.2 (solid line) and 0.1 (dashed line). 
The mass of the t quark was taken to be 170 GeV. 

FIG. 2. The mass of the bi squark is plotted as a function of the universal scalar mass mo for 
Rsgg = 0.3 (dot-dashed line), 0.2 (solid line) and 0.1 (dashed line). The mass of the t quark was 
taken to be 170 GeV. 
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